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Abstract. We present a Monte Carlo (MC) radiative transfer code for 
complex three dimensional dust distributions and include transiently 
heated PAH. The correctness of the code is confirmed by comparison 
with benchmark results. The method makes use of the parallelization 
capabilities of modern vectorized computing units like graphic cards. 
The computational speed grows linearly with the number of graphical 
processing units (GPU). On a conventional desktop PC, our code is 
up to a factor 100 faster when compared to other MC algorithms. As 
an example, we compute the dust emission of proto-planetary disks. 
We simulate how a mid-IR instrument mounted at a future 42m ELT 
will detect such disks. Two cases are distinguished: a homogeneous 
disk and a disk with an outward migrating planet, producing a gap 
and a spiral density wave. We find that the resulting mid-IR spectra 
of both disks are almost identical. However, they can be distinguished 
at those wavelengths by coronographic, dual-band imaging. Finally, 
the emission of PAHs exposed to different radiation fields is computed. 
We demonstrate that PAH emission depends not only on the strength 
but also strongly on the hardness of the radiation, a fact which has 
often been neglected in previous models. We find that hard photons 
(> 20 eV) easily dissociate all PAHs in the disks of T Tauri stars. To 
explain the low, but not negligible detection rate (< 10%) of PAHs in 
T Tau disks, we suggest that turbulent motions act as a possible path 
for PAH survival. 



1 Dust model 

In the dust model we consider: large (60A < a < 0.2 — 0.3/im) silicate (Draine 
2003) and amorphous carbon (Zubko et al. 2004) grains and small graphite (5 < 
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a < 80A) grains. For both, we apply a power law size distribution: n(a) oc a 3 5 
and absorption and scattering cross-sections are computed with Mie theory. In 
addition we include PAHs with 30 and 200 C atoms. For the absorption cross 
section (Draine (2010), Malloci (2010), Verstractc (2010)) we follow at photon 
energies between 1.7-15eV Malloci et al. (2007, their Fig. 4). At low frequency, we 
take a cross section cut-off given as an average of neutral (Schutte et al. 1993) and 
ionized (Salama et al. 1996) species. For hard photons, beyond the 17 eV band, 
we scale the photo-absorption cross section of PAH up to the keV region to similar 
sized graphite particles. By computing cross sections above 100 eV, we consider 
an approximation of kinetic energy losses (Dwek & Smith (1996) and apply it 
to all dust particles. With the advent of ISO and Spitzer more PAH emission 
features and more details of their band structures are detected (Tielens 2008). We 
consider 17 emission bands and take, for simplicity and as suggested by Boulanger 
et al. (1998) and Siebenmorgen et al. (1998), Lorentzian profiles. Parameters are 
listed in Siebenmorgen & Kriigel (2001) and calibrated using mid-IR spectra of 
starburst nuclei. For starburst galaxies, a SED library is computed (Siebenmorgen 
et al. 2007) and those models provide good fits to the SED of local galaxies and to 
PAH detected at high red shifts (z « 3, Efstathiou & Siebenmorgen 2009). In the 
model, we use dust abundances of [X]/[H] (ppm) of: 31 for [Si], 150 [amorphous 
C], 50 [graphite] and 30 [PAH], respectively; which is in agreement with cosmic 
abundance constrains (Asplund et al. 2009). We are in the process of upgrading 
the model to be consistent with the polarisation of the ISM (Voshchinnikov 2004). 

2 Monte Carlo method 

We compute the radiative transfer with a Monte Carlo (MC) technique, which 
allows to handle complicated geometries, by following the flight path of many 
random photons. The basic ideas of our method go back to Lucy (1999) and 
Bjorkmann & Wood (2001). They are described and extended to the treatment of 
transiently heated PAHs in Kriigel (2006). The space is partitioned into cubes and, 
where a finer grid is needed, the cubes are further divided into subcubes (cells). 
The star emits photon packages of constant energy e, but different frequencies. A 
package entering on its flight path a cell may be absorbed there or scattered. The 
probability for such an interaction is given by the extinction optical depth along 
the path within the cell. When the package is scattered, it only changes direction 
determined in a probabilistic manner by the phase function. When it is absorbed, 
a new package of the same energy, but usually different frequency v ncw is emitted 
from the spot of absorption. The emission is isotropic. Each absorption event 
raises the energy of the cell by e, and accordingly its temperature. 

We were able to increase the computational speed of the code by up to two 
orders of magnitude by parallelizing it, i.e. by calculating the flight paths of a 
hundred or more photons simultaneously. This allows the treatment of complex 
geometries which were so far prohibited because of their excessive computer de- 
mands. To achieve parallelization, we had to slightly modify the original code. 
The latter calculates the new frequency v ncvl of a package that is emitted after ab- 
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sorption using the cell temperature before absorption, but applies a correction for 
the new temperature due to Bjorkmann & Wood (2001). We now omit this correc- 
tion or in, mathematical terms, use Eqn. (11.67) instead of Eqn. (11.65) in Krugel 
(2006). As expected and borne out in tests, the correction may be neglected when 
the number of photons absorbed in a cell is not small. The vectorized MC code 
was developed during the PhD project of Heymann (2010). It was verified against 
benchmarks provided for stellar heated dust spheres by Ivezic et al. (1997), the 
ray tracing code by Krugel (2006) and axi-symmetrical disks by Pascucci et al. 
(2004). Generally, the computed SED agrees with the benchmarks to within a few 
percent. One example of such comparisons is shown in Fig. [T] for a dust sphere 
with PAHs. The factor by which parallelization speeds up the computation scales 
almost linearly with the number of graphical processing units. We point out that 
particular attention had to be given to the choice of the random number generator 
where we chose the Mersene Twister algorithm (Matsumoto & Nishimura, 1998). 
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Fig. 1. Comparison of a SED of a stellar heated dust sphere of constant density and 
visual extinction to the star of Av — 10 mag computed with a ray tracing method, as 
described in Krugel (2006), and the MC treatment of this work; both methods agree to 
within a few % . 

3 Detection of proto-planetary disk structures 

Hydrodynamical simulations of proto-planetary disks with an orbiting planet show 
particular disk features (Masset et al. 2007). We present a three dimensional 
application of the MC code where a disk heated by a solar type star has a dust 
density which follows the one as given by Pascucci et al., and has in addition a 
gap between 3 - 8AU and a low density spiral structure. The optical depth along 
the midplane from 0.2 to 75 AU is 10 mag. We simulate if such disk structures can 
be resolved at a distance of 50 pc with a mid-IR instrument (Brandl et al. 2010) 
mounted at a future 42m extreme large telescope (ELT); a project under study by 
ESO. The point spread function (PSF) of the ELT has resolution of 50 milliarcsec 
at 10/zm. In order to improve the contrast between star and disk, the instrument 
will provide coronographic and dual band imaging modes. We choose band passes 
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Fig. 2. Image at 11.3/im (left) together with flux profiles along the major axis (right). 
The 3a detection limit (dotted) and PSF (dash-dotted) is indicated. For the disk with 
(magenta) or without (black) PAH the gap and the spiral structure is well preserved. A 
homogeneous disk including PAHs is shown in dased for comparison. 



at 11.3 and 10/im. The 11.3/im emission is shown in Figj2] together with the 
flux profile along the major axis. Profiles of such a disk with and without PAHs 
and that of a homogeneous disk are computed. The 3(7 detection limit after lh 
integration is given assuming background limited performance of the instrument. 
The models predict that the detailed structure of such a proto-planetary disk can 
be well detected. 
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Fig. 3. The temperature distribution P(T) of a PAH with 100 C atoms exposed to mono- 
chromatic radiation with hv — 3.8, 13.6, 50, 100, 300 eV and 1 keV in a constant heating 
bath of flux F = 10 4 erg s _1 cm -2 (Siebenmorgen & Kriigel 2010). The shaded area 
marks temperatures above the sublimation temperature of graphite. 
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4 Destruction and survival of PAH in T Tauri disks 

Despite the fact that from the stellar heating one would expect to detect PAH in 
the disks of T Tauri, one find them rather seldom, in less than 10% of T Tauri 
(Geers et al., 2006). In order to explain this fact we present a simplified scheme to 
estimate the location from the primary heating source at which the PAH molecules 
become photo-stable (Siebenmorgen & Kriigel 2010). T Tauri stars have beside 
photospheric emission also a far ultraviolet (FUV), an extreme ultraviolet (EUV) 
and an X-ray component with a fractional luminosity of about 1%, 0.1% and 
0.025%, respectively. Such hard photons are very efficient in dissociating PAHs. 
The temperature distribution, P(T), of PAHs after hard photon absorption is 
shown in Fig. [3] It demonstrates that P(T) depends strongly on the hardness 
and spectral shape of the exciting radiation field, a fact which is often neglected in 
computations of the PAH emission. After photon absorption, a highly vibrationally 
excited PAH may relax through emission of IR photons or, if sufficiently excited, 
lose atoms (Omont (1986), and Tielens (2005) for a textbook description). We 
find that hard photons (EUV and X-ray) would destroy all PAHs in the disk of T 
Tauri stars; whereas soft photons with energies < 20 eV dissociate PAHs only up 
to short distances from the star (1-2 AU). As a possible path for PAH-survival 
turbulent vertical motions are suggested. They can replenish or remove PAHs 
from the reach of hard photons. In our treatment the presence of gas is considered 
which is ionized at the top of the disk and neutral at lower levels. A view of the 
scheme is shown as vertical cut along the midplane in Fig. |4] 
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Fig. 4. Of each radiation components (photosphere, FUV, EUV and X-rays) of the 
fiducial T Tauri star about 90 % is absorbed in what we call extinction layer. The 
optical depth from its bottom to the star is one and in the vertical direction equal to 
the grazing angle a s . The height of its lower boundary zo declines with the radius, but 
its geometrical thickness is rather constant (I « 0.5 H, for details see Siebenmorgen & 
Kriigel 2010). Vertical motions may either remove PAH from the extinction layer quickly 
enough for survival or replenish PAH from below. 
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